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Two-dimensional (2D) intrinsic ferromagnetic (FM) semiconductors are crucial to develop low¬ 
dimensional spintronic devices. Using density functional theory, we show that single-layer chromium 
trihalides (SLCTs) (CrX 3 ,X=F, Cl, Br and I) constitute a series of stable 2D intrinsic FM semicon¬ 
ductors. A free-standing SLOT can be easily exfoliated from the bulk crystal, due to a low cleavage 
energy and a high in-plane stiffness. Electronic structure calculations using the HSE06 functional 
indicate that both bulk and single-layer CrX 3 are half semiconductors with indirect gaps and their 
valence bands and conduction bands are fully spin-polarized in the same spin direction. The en¬ 
ergy gaps and absorption edges of CrBr 3 and CH 3 are found to be in the visible frequency range, 
which implies possible opt-electronic applications. Furthermore, SLCTs are found to possess a large 
magnetic moment of 3 pLB per formula unit and a sizable magnetic anisotropy energy. The magnetic 
exchange constants of SLCTs are then extracted using the Heisenberg spin Hamiltonian and the 
microscopic origins of the various exchange interactions are analyzed. A competition between a near 
90° FM superexchange and a direct antiferromagnetic (AFM) exchange results in a FM nearest- 
neighbour exchange interaction. The next and third nearest-neighbour exchange interactions are 
found to be FM and AFM respectively and this can be understood by the angle-dependent ex¬ 
tended Cr-X-X-Cr superexchange interaction. Moreover, the Curie temperatures of SLCTs are also 
predicted using Monte Carlo simulations and the values can further increase by applying a biaxial 
tensile strain. The unique combination of robust intrinsic ferromagnetism, half semiconductivity and 
large magnetic anisotropy energies renders the SLCTs as promising candidates for next-generation 
semiconductor spintronic applications. 


I. INTRODUCTION 

The successful isolation of graphene from graphite has 
triggered extensive experimental and theoretical studies 
on two-dimensional ( 2 D) materials. 1 Graphene has been 
found to possess many extraordinary properties including 
exceptionally high electron mobility ( ~ 10 5 cm 2 V _1 s _1 ) 
and Young’s modulus (~ 1 TPa). 2 Besides graphene, var¬ 
ious 2 D materials including hexagonal boron nitride, 3 
silicene 4 and transition metal dichalcogenides 5,6 such as 
M0S2 have also been identified as important candidates 
for next-generation electronic and optoelectronic devices 
due to their unique properties. 7 Some available 2 D ma¬ 
terials may also be combined to produce novel van der 
Waals (vdW) heterostructures, which may lead to even 
richer physics and device applications. 8 However, most 
of the available 2 D materials in the pristine form are 
intrinsically nonmagnetic, and this limits their applica¬ 
tions in spintronics. Although magnetic moments can 
be introduced by dopants, defects, edges, and coupling 
to FM substrates 9 , a long-range magnetic order is rarely 
observed experimentally in 2 D materials. It is thus desir¬ 
able to explore novel 2 D materials with a robust intrinsic 
ferromagnetic order. 

Layered transition metal compounds constitute a large 
family of materials with a wide range of electronic, op¬ 
tical, and mechanical properties. 6,10 The weak interlayer 
vdW bonding gives the possibility of their exfoliation into 
single-layer nanosheets. 10 Meanwhile, quantum confine¬ 


ment of electrons upon exfoliation is expected to lead 
to intriguing electronic, optical and magnetic proper¬ 
ties distinguished from their bulk. For example, bulk 
M0S2 is an indirect semiconductor with a band gap of 
1.2 eV, while single-layer M0S2 has a direct band gap of 
1.8 eV. 5,6 Due to the direct band gap in the visible fre¬ 
quency range, M0S2 is seen as an alternative to graphene 
and is considered as a promising candidate for next- 
generation electronic and optoelectronic devices. With 
unfilled d/f electrons in transition metal, the existence of 
magnetism is possible for transition metal compounds in 
both bulk and single-layer forms. Recently, there is also 
increasing interest for exploring novel single-layer tran¬ 
sition metal compound with long-range magnetic order. 
Various 2 D single-layer materials including K2CUF4 11 , 
doped-MnPS3 12 and CrSi(Ge)Te3 13 are predicted to be 
producible with an exfoliation process with the intrin¬ 
sic intralayer ferromagnetic order of the bulk. Moreover, 
due to the strong spin-orbit coupling of transition met¬ 
als, large magnetic anisotropy energies were also found in 
single-layer transition metal compounds 14 , and this holds 
promise for low-dimensional magneto-electronics applica¬ 
tions. 

In the present work, we focus on a series of layered 
semiconductors chromium trihalides, CrXs(X=F, Cl, Bl¬ 
and I) 15-17 . Although CrBr3 in particular was known as 
the first ferromagnetic semiconductor found in I 960 18 , 
the knowledge of CrX3 in general is still limited. It 
is known that CrX3(X=Cl, Br and I) crystalize in the 
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FIG. 1. (Color online) Crystal structure of ferromagnetic 
chromium trihalides CrX 3 in hexagonal (a) and rhombohe- 
dral (b) lattices. The red arrows represent the spin direction 
of Cr atoms. 

rhombohedral Bil 3 structure (space group R3 ) as shown 
in Fig. 1 at low temperature and can be transformed into 
the monoclinic AICI 3 structure (space group C2/m) upon 
warming. In both phases, the Cr 3+ ions are arranged in 
a honeycomb network, and coordinated by edge-sharing 
octahedra with six X - ions. The resulting halide-Cr- 
halide triple layers are stacked along the c-axis with a 
van der Waals gap. All three compounds are found to 
exhibit a strong intralayer ferromagnetic coupling, but 
the interlayer couplings seem to be distinct. CrCl 3 layers 
stacked antiferromagnetically with a Neel temperature 
(T^) of 16.8 K 19 , while CrBr 3 18 and Crl 3 17 layers are cou¬ 
pled ferromagnetically with Curie temperatures (Tq) of 
37 K and 61 K. In addition, the spins of the Cr atoms 
are found to lie in the CrX 3 plane, whereas the easy axes 
of the bromide and iodide are found to be along the c- 
axis. The combination of a weak interlayer coupling, an 
intrinsic magnetic order and semiconductivity of CrX 3 
provide an excellent platform to explore the existence 
of novel 2D magnetic semiconductors, which are needed 
by low-dimensional spintronics. Very recently, McGuire 
et al. 1 ' have re-examined the crystallographic and mag¬ 
netic properties of Crl 3 experimentally, and also inves¬ 
tigated both its bulk and single-layer properties using 
first-principles calculations. Their results indicate that 
Crl 3 layers can be easily exfoliated from bulk crystals 
and the ferromagnetic order is still retained in the single¬ 
layer form. However, detailed knowledge about the mag¬ 
netic interactions and magnetic anisotropy of chromium 
trihalides, which is crucial to understand the magnetic 
mechanism and potential applications, are still lacking. 

Here, we present a comprehensive theoretical inves¬ 
tigation on bulk chromium trihalides and their single¬ 
layer counterparts. The interlayer vdW interactions are 
included using the optB 88 -vdW functional, while the 
HSE06 functional is also used to obtain accurate elec¬ 


tronic structures and optical properties. The phonon 
dispersion calculations and ab initio molecular dynamics 
(MD) simulations are also performed to evaluate the sta¬ 
bility of single layer. Moreover, the magnetic exchange 
constants are extracted by comparing the total energy 
with the Heisenberg spin Hamiltonian. The Curie tem¬ 
perature is also predicted using Monte Carlo simulations. 
The microscopic origins of various exchange interactions 
involved are discussed in details. The obtained geome¬ 
tries, electronic structures and magnetic properties are 
then compared with available experimental results. Our 
calculations show that single-layer chromium trihalides 
are stable two-dimensional intrinsic ferromagnetic half¬ 
semiconductors with large magnetic anisotropy energies. 
These make them promising materials for future spin- 
tronic applications. 


II. COMPUTATIONAL DETAILS 

The present calculations have been performed with 
spin-polarization using the Vienna ab initio simulation 
package (VASP) code 20,21 within projector augmented- 
wave (PAW) method 22,23 . General gradient approxi¬ 
mations (GGA) in the Perdew-Burke-Ernzerhof(PBE) 
implementation 24 are applied as the exchange correlation 
functional. To better take into account the interlayer van 
der Waals forces of CrX 3 , we also perform calculations 
adopting non-local van der Waals density functional in 
form of optB 88 -vdW 25,26 , which has been evidenced to 
give a much improved descriptions for various systems in¬ 
cluding graphene/metal interfaces. 27 The range-separated 
hybrid functional HSE06 28 is used in calculating the elec¬ 
tronic structure and optical properties accurately. In 
the HSE06 method, the short-ranged exchange is con¬ 
structed by mixing 25 % exact nonlocal Hartree-Fock 
exchange and 75 % semi-local PBE exchange. While 
both the electron correlation and long-ranged exchange 
are still treated at the PBE level. Convergence tests have 
been performed carefully both for plane-wave cutoff en¬ 
ergy and fc-point sampling. A plane-wave basis set ex¬ 
panded in energy with a cutoff of 600 eV is used in the 
calculation. Jf-points sampling with a mesh of 8 x 8 x 8 k- 
points generated by the scheme of Monkhorst-Pack 29 is 
used for a primitive cell of bulk CrX 3 as shown in Fig. 1- 
(b). And a 8x8x1 mesh is used for a single-layer CrX 3 
(Fig. 3). All the lattice constants and atomic coordinates 
are optimized until the maximum force on all atoms is 
less than 5 x 10 -3 eV/A. A large vacuum space of at 
least 18 A thick is included in the supercell to avoid 
interaction between images. The elastic constants are 
obtained by an energy-strain approach, in which defor¬ 
mations applied to the unit cell and the corresponding 
energy-strain relationship can be found in Ref. 30 . The 
phonon calculations have been performed using the finite 
displacement approach, as implemented in the Phonopy 
code 31 , in which a 3x3x1 superccll and a displacement 
of 0.01 A from the equilibrium atomic positions are em- 
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ployed. 

To extract the exchange interactions parameters be¬ 
tween spins, we fit the total energies from DFT calcu¬ 
lations for various spin configurations to the Heisenberg 
spin Hamiltonian 

H=-J 1 J2s i ■ s 0 Sj -Js Sj 

<*i) <(b» <(<b’») 

(!) 

defined on a honeycomb lattice. Here, |/>|=3/2. Using 
these DFT-derived magnetic exchange parameters, the 
Curie temperature is then estimated using Metropolis 
Monte Carlo (MC) simulations, in which a 100 xlOO x 1 
2D honeycomb lattice with periodic boundary conditions 
is used. For each temperature studied, the MC simula¬ 
tion involves for 10 5 MC steps per site to attain thermal 
equilibrium. The temperature-dependent mean magnetic 
moment and magnetic susceptibility are then obtained. 
The temperature at which the mean magnetic moment 
drastically drops to nearly zero and the magnetic suscep¬ 
tibility peaks sharply is identified as the Curie tempera¬ 
ture. 


III. RESULTS AND DISCUSSION 
A. Structure, Cleavage and Stability 

First, we discuss the structural properties of bulk 
chromium trihalides. CrX 3 (X=Cl,Br and I) crystallize 
in the R3 phase at low temperature, and the correspond¬ 
ing hexagonal and rhombohedral lattices are shown in 
Fig. 1. It should be pointed out that CrF 3 is scarcely 
discussed in the literature. An experiment 32 nevertheless 
indicates that it takes the R3c phase. From experiments, 
the interlayer coupling is ferromagnetic for CrBr 3 and 
Crl 3 , but antiferromagnetic for CrCl 3 . Both PBE and 
optB88-vdW calculations however indicate that the ferro¬ 
magnetic state is more stable than the non-magnetic and 
anti-ferromagnetic phases in bulk for all four chromium 
trihalides. Similar results were also reported in early 
calculations. 33 Since the Cr atoms in the CrX 3 plane 
couple ferromagnetically in both ferromagnetic and anti¬ 
ferromagnetic bulk CrX 3 , the inconsistency of the inter¬ 
layer magnetic state between experiment and theory for 
CrCl 3 should not effect the following calculation about 
single-layer chromium trihalides. We thus only discuss 
the results of ferromagnetic CrX 3 . 

The lattice constants of the ferromagnetic CrX 3 cal¬ 
culated using the PBE and optB88-vdW functionals are 
given in Table. I. Experimental results available from 
the literature are also listed for comparison. As shown 
in the table, both lattice parameters a and c increase 
with the element number of X. This can be attributed 
to the increasing atomic radius of X and the weakening 
reactivity between the X and the Cr atoms. It is clearly 
seen that the PBE functional significantly overestimates 
the interlayer distance c by about 10%. When the van 


TABLE I. The calculated lattice parameters of CrX 3 using 
PBE and optB88-vdW functionals. Experimental results are 
also listed for comparison. For CrF 3 , the lattice parameters 
of R3c phase are used. 




PBE 

optB88-vdW 

experiment 

CrF 3 

a (A) 

5.185 

5.163 

4.992“ 


c(A) 

14.995 

13.351 

13.215 

CrCl 3 

a (A) 

6.053 

5.987 

5.942'’ 


c(A) 

19.422 

17.153 

17.333 

CrBi '3 

a (A) 

6.435 

6.350 

6.260“ 


c(A) 

20.836 

18.259 

18.200 

Crl 3 

a (A) 

7.006 

6.905 

6.867 d 


c(A) 

22.350 

19.805 

19.807 


“Ref. 32 
6 Ref. 15 
“Ref. 16 
d Ref. 17 


der Waals forces are treated using the optB88-vdW func¬ 
tional, a good agreement between theory and experiment 
is achieved. This demonstrates that vdW forces play a 
very important role in interlayer binding and the optB88- 
vdW functional provides accurate results. It will thus 
be used in calculating the cleavage energy of single-layer 
chromium trihalides. 

We now focus on the possibility of the exfoliation of 
chromium trihalide layers from the bulk and the stability 
of free-standing single layers. To obtain a free-standing 
membrane, the cleavage energy needed to be overcome 
in the exfoliation process should be small. In addition, 
a high in-plane stiffness is also needed to avoid curling 
upon exfoliation. 11 . A large gap d between two layers 
representing a fracture in the bulk as shown in Fig. 2-(a) 
is introduced to simulate the exfoliation procedure. The 
cleavage energies as functions of d calculated using the 
optB88-vdW functional are shown in Fig. 2-(b). The 
energy relative to the un-cleaved equilibrium state in¬ 
creases with the separation d, and gradually converges 
to the cleavage energy. It can be seen that the cleav¬ 
age energies for the four CrX 3 materials are very simi¬ 
lar and ranges only from 0.28 to 0.30 J/m 2 . These val¬ 
ues are comparable with other similar systems: CrSiTe 3 
(0.35 J/m 2 ), CrGeTe 3 (0.38 J/m 2 ) 13 , and MnPSe 3 (0.24 
J/m 2 ) 12 . They are also smaller than the experimentally 
estimated cleavage energy for graphite (0.36 J/m 2 ). This 
indicates that the exfoliation of single layers from bulk 
chromium trihalides is feasible experimentally. The the¬ 
oretical cleavage strength given by the maximum deriva¬ 
tive of the cleavage energy is 1.495, 1.384 ,1.303 and 1.177 
GPa for CrF 3 ,CrCl 3 , CrBr 3 and Crl 3 , respectively. 

In addition, the in-plane elastic constants of single¬ 
layer chromium trihalides are also calculated using the 
PBE functional. The primitive cell of a 2D CrX 3 layer 
used in the calculations is shown in Fig. 3. The cal¬ 
culated in-plane lattice constants as shown in Table. II 
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FIG. 2. (Color online) Supercell model (a) of CrX 3 used to 
simulate the exfoliation procedure and the cleavage energy (b) 
calculated using the optB88-vdW functional as a function of 
the separation (d) between two fractured parts, do represents 
the equilibrium interlayer distance of chromium trihalides. 




FIG. 3. (Color online) The atomic structure and unit cell 
of ferromagnetic single-layer chromium trihalides upon exfo¬ 
liation used in the electronic and phonon calculation. The 
transparent blue and red atoms represent the X atoms in top 
and bottom layer. 


are very similar to the bulk values. Moreover, cn, C 12 , 
Young modulus Y 2 d and Poisson ratio v decrease mono- 
tonically from F to I, which is also consistent with the 
binding strength between Cr 3+ and X~. We further es¬ 
timate the gravity induced out-of-plane deformation us¬ 
ing the formula h/L « pgL/Y^D , 34 where the length of 
the flake L is set to 100 pm, and p is the density of 2D 
chromium trihalides. The resulting values of h/L for 2D 
CrX 3 (X=F,Cl,Br and I) are 3.405x 10~ 4 ,3.629x 10~ 4 , 
4.488xl0 -4 , and 5.178xl0 -4 . For graphene, h/L is 
about 10 -4 for flakes of this size. This suggests that 
free standing single-layer chromium trihalides are able to 
withstand their own weight and not curl easily. The low 
cleavage energy and strong in-plane stiffness promise that 
the 2D chromium trihalides can be produced by a simple 


FIG. 4. (Color online) Phonon dispersion of single-layer 
chromium trihalides calculated using PBE functional. 


exfoliation process. 

To further confirm the stability of free-standing single¬ 
layer chromium trihalides, we have calculated their 
phonon dispersion. As shown in Fig. 4, the imaginary 
frequency is found to be absent in the whole Brillouin 
Zone for all 2D chromium trihalides. This suggests that 
these single layer materials are dynamically stable and 
can exist as freestanding 2D crystals. In addition, the 
phonon frequency is found to be softening monotonically 
from F to I due to the increasing atomic weight of X. 
Furthermore, we have performed ab initio molecular dy¬ 
namics simulations at 300 K in canonical ensemble using 
Nose heat bath scheme to examine the thermal stability 
of the single-layer chromium trihalides. In these calcula¬ 
tions, a large supercell of (4x4) is used to minimize the 
constraint of periodic boundary condition. Our results 
indicate that the atomic configurations of the single-layer 
chromium trihalides including the honeycomb network of 
the Cr 3+ ions and octahedra of CrXg remain nearly intact 
after heated for more than 9 ps with a time step of 3 fs. 
The variations of the total potential energy with respect 
to simulation time are plotted in Fig. 5. As shown in the 
figure, the total potential energies remain almost invari¬ 
ant during the simulation, which suggests that the single¬ 
layer chromium trihalides are thermally stable at room 
temperature. Clearly, the low cleavage energy, dynamics 
stability and thermal stability of single-layer CrX 3 pre¬ 
dicted here suggest that the free-standing ferromagnetic 
chromium trihalides single layers can be realized experi¬ 
mentally even at room temperature. 
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TABLE II. The structural, elastic properties, cleavage strength(cr) magnetic anisotropy energy (MAE) and easy axis (EA) of 
single-layer chromium trihalides calculated using PBE functional. 



°(A) 

cii (Nm 1 ) 

C 12 (Nm -1 ) 

YjofNin- 1 ) 

V 

P(Kg m 2 ) 

h/L 

g (GPa) 

MAE (fj,eV/Cr) 

EA 

CrF 3 

5.189 

49.947 

23.862 

38.547 

0.478 

1.553X 10 -6 

3.405X10" 4 

1.495 

119.0 

c 

CrCl 3 

6.051 

37.308 

11.069 

34.024 

0.297 

1.659X 10 -6 

3.629X 10 -4 

1.384 

31.5 

c 

CrBr 3 

6.433 

31.741 

8.811 

29.295 

0.278 

2.703X10 -6 

4.488X 10 -4 

1.303 

185.5 

c 

Crl 3 

7.008 

25.474 

6.443 

23.845 

0.253 

3.378X10 -6 

5.178X 10 -4 

1.177 

685.5 

c 



(a):CrF 3 



Time (fs) 
(c):CrBr 3 



2000 4000 6000 8000 


Time (fs) 

(b):CrCl 3 



2000 4000 6000 8000 
Time (fs) 

(d):CrI 3 


FIG. 5. (Color online) Variations of the total potential energy 
of single-layer chromium trihalides with respect to simulation 
time during ab initio molecular dynamics simulations. 


B. Electronic Structures and Optical properties 

We further investigate the electronic structures of 
single-layer chromium trihalides. Since PBE functional 
usually underestimates the energy band gap, hybrid func¬ 
tional in form of HSE06 functional is used to get accurate 
electronic structures. As shown in Fig. 6a, all monolayer 
chromium trihalides are indirect semiconductors. The 
band gap decreases with element number of X and the 
corresponding values are 4.68, 3.44, 2.54, and 1.53 eV 
for CrX 3 (X=F,Cl,Br and I) respectively. Moreover, we 
also find from Fig. 7a that the valence band and conduc¬ 
tion band edges around the Fermi level are fully spin- 
polarized and exclusively contributed by same spin com¬ 
ponent, which shows a typical half-semiconductor charac¬ 
ter. The difference of the band edge energy between the 
two spin components (denoted by A cbm = E%°™ n - E^ m 
and A vbm = E^™ n — E^ m for the conduction band min¬ 
imum and the valance band maximum respectively) are 
also reported in Tab. III. In particular, A cbm reaches 


1.609 eV in the case of Crl 3 , which implies possible ap¬ 
plications in spin-polarized carrier injection and detec¬ 
tion. More interestingly, occupied Cr-3d orbitals are only 
found in the spin-up direction and the spin-down Cr-3d 
states are fully unoccupied. These results can be under¬ 
stood using a crystal field theory. Since the Cr 3+ ions 
are located at a octahedra environment coordinated by 
six X - ions, the Cr-3d orbitals are spitted into three low- 
lying t 2 g and two high-lying e g orbitals. According to 
Hund’s rule and the Pauli exclusion principle, Cr 3+ will 
take a high-spin t 3 g e 3 electronic configuration, which 
implies that the occupied 3d orbitals must all be in one 
of the spin directions. Moreover, we can also find in¬ 
teresting element- and spin-dependent gap character in 
CrX 3 . As shown in Fig. 7a, HSE06 calculations indi¬ 
cate that the conduction band is dominated by Cr-3d 
states weakly hybridized with the X-3p states for all four 
compounds in both spin directions. And the states of 
the valence band in spin-down direction are of almost 
pure X-3p character, which is also element-independent. 
For the upper part of the valence band, the states are 
contributed by the X-3p states with a mixture of Cr-3d 
states in CrX 3 (X=Cl,Br and I) and the contribution of 
X-3p increases with the element number of X. However, 
the Cr-3d states of CrF 3 still dominate the valence band 
near fermi level. Hence the band gaps in spin-down di¬ 
rection predicted by HSE06 are of charge-transfer-type 
but shows a element-dependent behaviour in spin-up di¬ 
rection. We can find that the gap is Mott-Hubbard type 
for CrF 3 and will gradually change from a mixture of a 
Charge-Transfer and a Mott-Hubbard type for CrCl 3 and 
CrBr 3 to Charge-Transfer type in case of Crl 3 . Results 
based on the PBE functional (Fig. 6b and Fig. 7b) also 
implies indirect semiconductors but with a much smaller 
band gap. However, it should be noticed that HSE06 and 
PBE functionals give different hybridization behaviours. 
Compared with HSE06 results, the crystal field splitting 
and exchange splitting of the Cr-3d states are strongly de¬ 
ceased in PBE functional (Fig. 7b). And the Cr-3d states 
(i.e. t 2 g states) are shifted to higher energies closed to 
Fermi level compared with the X-3p states and predomi¬ 
nated the valence band maximum, resulting in a strongly 
increased hybridization between X-3p and Cr-3d states. 
PBE functional thus tends to give a Mott-Hubbard type 
band gap in spin-up direction whereas the gap character 
in spin-down direction is still same to HSE06. Further 
experiments should be conducted to clarify this interest¬ 
ing spin-dependent electronic structure of CrX 3 . 
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(a) (b) 

FIG. 6. (Color online) Band structures of single-layer chromium trihalides calculated using HSE06 (a) and PBE (b) functionals. 
The red(blue) lines represent the band structure in the spin-up (spin-down) direction. 



- Total - X-p - Cr-d 






CrBr3 Crl3 


(a) (b) 

FIG. 7. (Color online) Total and partial density of states of single-layer chromium trihalides calculated using HSE06 (a) and 
PBE(b) functionals. 


We have also calculated the optical properties of CrX 3 . 
As shown in Fig. 8a, due to hexagonal symmetry of the 
CrX 3 lattice, the dielectric functions along the x and y 
axis are identical but different from that along the z axis, 
implying anisotropic optical properties. We can also find 
that the main peak of e 2 in CrX 3 can be traced back 


to the p-d excitation from X-3p states of valence band 
to Cr-3d states of lowest conduction band. For CrF 3 , 
the satellite states around 6 eV are contributed by d-d 
excitation from occupied 3d states around -1 eV and un¬ 
occupied states around 5 eV as shown in Fig. 7a. Our 
PBE calculations (Fig. 8b) give similar results but the ab- 
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TABLE III. The characteristic parameter in unit of eV of energy band of single layer and bulk chromium trihalides calculated 
using PBE and HSE06 functionals. 






PBE 








bulk 




single layer 



^cbm 

^vbm 

Gap up 

G&Pdown 

^cbm 

^vbm 

Gapnp 

G&P down 

CrF 3 

0.816 

-3.102 

1.609 

5.527 

0.818 

-3.146 

1.614 

5.578 

CrCl 3 

0.474 

-1.324 

1.454 

3.252 

0.451 

-1.309 

1.521 

3.281 

CrBr 3 

0.62 

-0.721 

1.182 

2.523 

0.624 

-0.776 

1.331 

2.731 

Crl 3 

0.754 

-0.231 

0.913 

1.898 

0.686 

-0.280 

1.143 

2.109 





HSE06 








bulk 




single layer 



^cbm 

^vbm 

G&piip 

Gap down 

^cbm 

^vbm 

Gap^p 

Gap down 

CrF 3 

1.075 

-1.862 

4.728 

7.665 

1.14 

-2.071 

4.680 

7.891 

CrCl 3 

0.987 

-0.267 

3.305 

4.559 

0.992 

-0.343 

3.441 

4.776 

CrBr 3 

1.316 

-0.167 

2.402 

3.885 

1.322 

-0.220 

2.538 

4.080 

Crl 3 

1.639 

-0.204 

1.286 

3.129 

1.609 

-0.243 

1.525 

3.377 



(a) (b) 

FIG. 8. (Color online) Imaginary part of the dielectric function t 2 of single-layer chromium trihalides calculated using HSE06 
(a) and PBE (b) functionals. 


sorption edges at lower energies than those found using 
HSE06. In addition, the stronger satellite states induced 
by d-d excitation are found to exist in CrX 3 , which can be 
understood by the fact that the valence band maximum 
is predominated by Cr-3d states in PBE calculation. It 
should be noticed that the main absorption peaks for 
CrBi '3 and Crl 3 fall in the visible range, although they 
are indirect semiconductors. This suggests that CrBr 3 
and Crl 3 are possible candidate for optoelectronic appli¬ 
cation. 

To reveal the change of electronic and optical proper¬ 
ties of CrX 3 due to quantum confinement, we have also 


calculated the electronic structure and optical proper¬ 
ties of bulk CrX 3 using the HSE06 and PBE functionals. 
Our results indicate that the band structure (Fig. 9) and 
density of states (Fig. 10) of bulk CrX 3 are very similar 
to their single layer counterparts. This is except that 
the band gap is slightly smaller, which is also found in 
most 2D transition metal compounds. Experimentally, 
the band gaps of bulk CrBi -3 and Crl 3 are 2.1 eV 35 and 
1.2 eV 36 respectively, which are in good agreement with 
2.40 and 1.28 eV given by our HSE06 calculations. We 
also compare our calculated imaginary part of the di¬ 
electric function with available experimental results for 




















CrBi '3 and CrCl 3 37 . As shown in Fig. lib, PBE under¬ 
estimates the band gap significantly, and this leads to an 
absorption edge at a lower energy than the experimental 
value. In addition, PBE also gives a satellite peak which 
seems to be invisible in experiment. The HSE06 calcu¬ 
lations seem to give more accurately predictions of the 
optical properties and a good agreement between theory 
and experiment is achieved in cases of CrBr 3 and CrCl 3 . 
This also verifies that our HSE06 calculations are reliable 
and are hence expected to give accurate band gaps and 
optical properties for single-layer chromium trihalides. 


C. Magnetic Properties 

We now evaluate the magnetic ground state of single- 
layer CrX 3 . Four possible different magnetic configura¬ 
tions as shown in Fig. 12 are considered 38 . Our calcu¬ 
lations indicate that the most stable magnetic state is 
the FM state. This means that intra-ferromagnetism in 
the bulk form are retained upon exfoliation. The mag¬ 
netic moment of a single unit cell is equal to 3 ac¬ 
curately and the local magnetic moments of Cr in CrX 3 
calculated using the PBE (HSE06) functional are 2.871 
(2.900), 2.937 (3.021), 3.012 (3.127), and 3.103 (3.305) 
Hb ■ This suggests that single layer CrX 3 are robust in¬ 
trinsic ferromagnetic materials with large magnetic mo¬ 
ments. 

Magnetic anisotropy energy (MAE) is another tech¬ 
nologically important parameter for magneto-electronic 
applications. Ferromagnetic materials with small MAE 
will result in super-paramagnetic rather than ferromag¬ 
netic behavior. We have thus calculated the MAE of 
single-layer chromium trihalides using the PBE func¬ 
tional which includes the spin-orbit coupling (SOC). The 
results as shown in Table. II indicate that easy axis is 
along c-direction for all CrX 3 compounds and isotropy 
holds in the basal plane. The corresponding MAE is 120, 
32,186, and 686 /xeV per Cr atom for CrF 3 , CrCl 3 , CrBr 3 
and Crl 3 , respectively. Bulk &I 3 was also reported to 
have the strongest magnetic anisotropy among chromium 
trihalides , 17 which is also in accordance with our predic¬ 
tion. It should be pointed out that the spins of the Cr 
atoms are found to lie in the basal plane in bulk CrCl 3 , 
whereas bulk CrBr 3 and Crl 3 have an easy axis along 
the c-direction. To confirm the easy axis of single-layer 
CrCl 3 , we have also calculated magnetic anisotropy en¬ 
ergy of single-layer &CI 3 using HSE06 functional. Our 
HSE06 results indicate that the easy axis of single-layer 
CrCl 3 is c-axis with a magnetic anisotropy energy 13.5 
/mV per Cr atom, which is in consistent with our PBE 
calculations. The transition from an easy-plane to an 
easy-axis for &CI 3 upon exfoliation could be under¬ 
stood by interplay between crystal field splitting and 
spin-orbital coupling, which is beyond the scope of the 
present work. We also find that the MAE of CrBi -3 and 
Crl 3 are much larger than those of commonly used ferro¬ 
magnetic materials such as bulk Fe (1.4 /reV per atom), 


and Ni (2.7 /reV per atom) 39 , and are also comparable 
to those of Ca 2+ doped La(Mn,Zn)AsO (509 /xeV per 
atom) 14 . The large MAE of single-layer CrBi -3 and Crl 3 
predicted here makes them promising candidates for low¬ 
dimensional magneto-electronic applications. 

In order to calculate the magnetic exchange interac¬ 
tions, we considered four different magnetic configura¬ 
tions as shown in Fig. 12. Using Eq. (1), the magnetic 
energy of two CrX 3 formulas can be explicitly expressed 
as: 

EpM/Neel = E 0 — (±3Ji + 6 J 2 ± 3J 3 )|S ? | 2 (2) 

Ezigzag/stripy = E 0 - (±Ji - 2 J 2 =F 3J 3 )|£| 2 (3) 

We then obtain exchange couplings by comparing with 
DFT results and also the Curie temperature using MC 
simulations. As listed in Table IV, both the nearest 
neighbour (NN) exchange interaction Ji and the next 
nearest neighbour ( NNN ) interaction J 2 are ferromag¬ 
netic, but the third nearest neighbour (3 NN) interac¬ 
tion J 3 is anti-ferromagnetic. Both the magnitudes of 
the magnetic interaction and the Curie temperature in¬ 
crease from F to I. To evaluate the importance of J 2 and 
J 3 , we also calculate Jnn only considering the NN in¬ 
teractions. Our calculations indicate that Jnn is slightly 
smaller than Ji. However, it should be noted that the 
magnitude of J 2 is one order larger than J 3 and the 
weight coefficient of J 2 in the Hamiltonian is twice that 
of Ji in the Ferromagnetic state. In particular, for CH 3 , 
we find that the ratio of the contribution of J 2 and J\ 
in the Hamiltonian as evaluated by 2 is about 0.45, 
which implies the importance of J 2 in the accurate de¬ 
termination of the Curie temperature. Using JnnJ^c of 
chromium trihalides are estimated to 36 ,36 ,51 and 56 K 
respectively. The corresponding experimental Curie tem¬ 
perature of bulk CrBr 3 18 and Crl 3 17 are 37 K and 61 K, 
respectively. When the next and third nearest neighbour 
exchange interactions are included, the T c increases to 
41, 49, 73 and 95 K respectively. To evaluate the ac¬ 
curacy of the PBE calculations for magnetic exchange 
interactions, we also calculate the exchange interaction 
only considering the NN interaction using HSE06, which 
is also listed in the Table. IV for comparison. The results 
given by HSE06 turn out to be very close to the PBE ex¬ 
cept the CrF 3 , which indicates the validity of the present 
PBE calculations. 

Now, we turn our attention to the microscopic origin 
of the exchange interactions in the single-layer chromium 
trihalides. The NN interaction Ji can be well understood 
from the competition between the direct exchange of the 
Cr-Cr sites and the superexchange mediated through the 
X ions. 38,40 . As shown above, the electronic configura¬ 
tion of Cr 3+ is t 2 g e° g . The corresponding direct exchange 
will be weekly antiferromagnetic 41 . Taking CH 3 as an 
example, the Cr-X-Cr angle (95.4°) is close to 90° in 
superexchange path shown in Fig. 13-(a), which is ex¬ 
pected to be ferromagnetic. Due to the large distance 
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(a) (b) 

FIG. 9. (Color online) Band structure of bulk chromium trihalides calculated using HSE06 (a) and PBE (b) functionals. The 
red(blue) lines represent the band structure in the spin-up (spin-down) direction. 


- Total - X-p - Cr-d 



- Total - X-p - Cr-d 






CrBr3 Crl3 


(a) (b) 

FIG. 10. (Color online) Total and partial density of states of bulk chromium trihalides calculated using HSE06 (a) and PBE 


(b) functionals. 

(4.05 A) between the two Cr ions, which reduces direct 
antiferromagnetic exchange, the FM superexchange will 
dominate the total exchange interaction. The underlying 
mechanism of J 2 and J 3 is more subtle. The cations are 
separated by two anions in the possible superexchange 
path, forming an extended cation-anion-anion-cation in¬ 


teraction path. Goodenough argued 41 that the rules 
for an extended cation-anion-anion-cation path are the 
same as those for cation-anion-cation interactions. The 
extended cation-anion-anion-cation interaction can thus 
be seen as the sum of exchange interactions from two 
cation-anion-anion ” superexchange path”. In addition, 
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- xx(yy) - zz - Exp. 


- xx(yy) - zz - Exp. 







(a) 


(b) 


FIG. 11. (Color online) Imaginary part of the dielectric function of bulk chromium trihalides calculated using HSE06 (a) and 
PBE (b) functionals. Available experimental data from Ref. 37 are also shown for comparison. 




FIG. 12. (Color online) The possible magnetic configuration 
of single layer chromium trihalides CrX 3 : FM (a), AF-Neel 
(b), AF-zigzag (c), and AF-stripy ordered (d). The red arrows 
represent the spin direction of Cr atoms. The crystal cells 
used in the calculation are also shown in Figure. 


TABLE IV. The exchange coupling and Curie temperature of 
single-layer chromium trihalides. 



tHSE06 

J NN 

meV 

Jnn 

meV 

Ji 

meV 

J 2 

meV 

J 3 

meV 

rpNN 

K 

T c 

K 

CrF 3 

0.8688 

1.7904 

1.8028 

0.0672 

-0.0128 

36 

41 

CrCl 3 

1.8728 

1.7908 

1.9232 

0.2264 

-0.1328 

36 

49 

CrBrs 

2.4216 

2.4524 

2.5988 

0.3800 

-0.1464 

51 

73 

Crl 3 

3.3216 

2.7096 

2.8624 

0.6384 

-0.1532 

56 

95 



according to the Goodenough-Kanamori-Anderson rule 41 
which has been used widely to understand the signs and 
strengths of superexchange interactions, a 180° superex¬ 
change is antiferromagnetic while a 90° superexchange 


FIG. 13. (Color online) Possible paths for the mediation of 
NN (a), NNN (b) and 3NN (c) magnetic exchange interac¬ 
tions in ferromagnetic single-layer Crl 3 . (Anti-)Parallel ar¬ 
rangement of the Cr 3+ ions represent the (anti-)ferromagnetic 
exchange interaction. 
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• Jnn ■—■ J\ □ Jl ®—® h ■—■ T c 



(c):CrBr3 (d):Crl3 


FIG. 14. (Color online) The strain dependence of exchange 
interactions and Curie temperature of single-layer chromium 
trihalides. 

interaction is ferromagnetic. For an intermediate angle 
between 90° and 180° , FM and AFM will compete with 
each other and a crossover angle from FM to AFM su¬ 
perexchange should exist. Subramanian et a /. 42 have in¬ 
vestigated the FM-AFM transformation by controlling 
the cation-anion-cation angle in cuprate perovskite and 
find a crossover angle about 127±0.6°. Thus, we can 
reasonably assume that a large cation-anion-anion angle 
( >130° ) contributes a AFM interaction while an angle 
near 90° presents a FM contribution. As suggested by 
Sivadas et al. 3S , there are many possible exchange paths 
for Ji and J 3 . Since the shortest distance between the 
X - pairs in different layers is larger than that between 
in-plane pairs, the most possible exchange path for J 2 
and J 3 should involve two X anions on the same plane. 
As shown in Fig. 13-(c), both cation-anion-anion angle 
for J 3 is 132.7° . Thus, J 3 should be AFM due to both 
AFM contributions from the two cation-anion-anion su¬ 
perexchanges. However, for the NNN exchange interac¬ 
tion in Fig. 13-(b), the two cation-anion-anion angles are 
139.8° and 97.0°, which contributes AFM and FM su¬ 
perexchange interactions, respectively. The sign of J 2 
will depend on the competition of the FM and AFM su¬ 
perexchange interactions. In CrX 3 , FM dominates over 
AFM and the resulting extended superexchange is FM. 
In contrast, for other systems such as M 11 PS 3 (Se 3 ) and 
CrSi(Ge)Te 3 , the total NNN exchange interaction with 
the extended superexchange path is AFM . 38 

Furthermore, we also investigate the magnetic inter¬ 
actions and Curie temperature of &X 3 under a biax¬ 
ial strain. Results obtained based on interactions up to 
nearest-neighbours or third-nearest neighbours are shown 
in Fig. 14. We find that the NN and NNN exchange in¬ 


teractions and the corresponding Curie temperature in¬ 
crease with strain, whereas the magnitude of 3NN inter¬ 
action decrease. Yoshida et al. 43 have investigated ex¬ 
perimentally the pressure effect on the Curie tempera¬ 
ture of CrBr 3 and found that T c decreases linearly with 
pressure, which agrees with our calculations. The Curie 
temperature can be tuned gradually with strain and 130 
K can be reached in the case of &I 3 using a 10% strain. 
Similar results are also reported in other systems such as 
CrSiTe 3. 13,38,40 In addition, we also find that a compres¬ 
sive strain can also induce a FM to AFM transition. The 
AF-Neel phase becomes the most stable phase at a strain 
of 6 %, 4%, 6 %, 8 % for CrF 3 , &CI 3 , CrBr 3 and CU 3 , re¬ 
spectively. Such a large compressive strain required to 
destroy the ferromagnetic state implies the robustness of 
the intrinsic ferromagnetism in single-layer &X 3 . 


IV. CONCLUSIONS 


In summary, we have investigated the geometry, elec¬ 
tronic, optical and magnetic properties of chromium tri¬ 
halides in single-layer and bulk forms. Our results in¬ 
dicate that the optB 88 -vdW functional can capture the 
interlayer vdW interaction of bulk &X 3 quite well while 
HSE06 gives consistent electronic structures and opti¬ 
cal properties with available experiment. Single-layer 
chromium trihalides are found to have very low cleav¬ 
age energies and high in-plane stiffness, which suggests 
that free-standing single-layer chromium trihalides can 
be easily exfoliated from the bulk crystals. The stability 
of single-layer chromium trihalides is further confirmed 
using the phonon dispersions calculation and ab initio 
molecular dynamics simulations. Furthermore, our re¬ 
sults also suggest that these 2D crystals exhibit robust 
long-range ferromagnetic order with magnetic moments 
of 3 hb per formula unit and large magnetic anisotropy 
energies. The ferromagnetic nearest-neighbour exchange 
interactions can be attributed to the competition be¬ 
tween direct antiferromagnetic exchange interactions and 
near 90° ferromagnetic superexchange interactions. The 
next-nearest-neighbour FM exchange interactions can be 
understood by the angle-dependent extended superex¬ 
change Cr-X-X-Cr interaction. Monte Carlo simulations 
based on the classical Heisenberg model predicts Curie 
temperatures of chromium trihalides up to 90 K, and can 
be further increased by applying biaxial tensile strains. 
Electronic structure calculations indicate that both bulk 
and 2D chromium trihalides are half semiconductors with 
indirect band gaps and the absorption edges of CrBi -3 and 
Crl 3 lie in the visible range, which suggests that they are 
promising candidates for future low-dimensional semi¬ 
conductor spintronic, magneto-electronic and magneto¬ 
optic applications. 
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